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Abstract 

A mathematical theory is developed to perform the calculations necessary to determine the wave drag 
for slender bodies of non-circular cross section. The derivations presented in this report are based on 
extensions to supersonic linearized small perturbation theory. A numerical scheme is presented utilizing 
Fourier decomposition to compute the pressure coefficient on and about a slender body of arbitrary cross 
section. 


Introduction 

Supersonic linearized small perturbation theory has been used to calculate the wave drag for simple aircraft 
configurations of slender fuselages and thin wings since the 1950s. Complex aircraft configurations have been 
analyzed by linearized small perturbation theory since the 1970s. The interference effects between aircraft 
fuselage, wing, and engine nacelle found in complex configurations can be simulated by the superposition 
of pressure fields from and on interfering bodies. Only axisymmetric slender body theory has been used 
to model the interference of aircraft fuselages and engine nacelles in the past. However, current designs 
for supersonic transport aircraft incorporate engine nacelle configurations which are not axisymmetric. In 
order to analyze the engine/airframe interference effects for these current nacelle designs other methods have 
been used. In particular, Euler Computational Fluid Dynamic techniques have been used to compute the 
interference wave drags for complex aircraft configurations. Using CFD techniques greatly complicates the 
analysis set up and run times, thus limiting the number of possible configurations which can be analyzed. 
The motivation for this work is to find an extended solution technique whereby linearized small perturbation 
theory can be used to calculate the pressure field about slender, non-axisymmetric bodies. This would extend 
the modeling capacities of existing aerodynamic configuration analysis programs based on linearized small 
perturbation methods. 

Governing Differential Equation and General Solution 

The following are the basic relations and theory developed in the non-axisymmetric wave drag analysis. The 
problem is formulated in cylindrical coordinates, as is typical for slender body theories. Figure 1 shows 
a slender body whose cross sections is not axisymmetric, and the coordinate system used in the following 
analyses. The x axis is aligned with the free stream flow direction. The governing partial differential equation 
for the perturbation velocity potential is: 


4 . + l d A+L&± 

l 1 M ^ dx 2 + dr 2 + rflr + r 2 <90 2 


( 1 ) 


Implicit in the derivation of Equation 1 are the assumptions that the gas flow is; supersonic (M > 1), ideal, 
irrotational, and inviscid. Ward (Reference 1) gives a general solution of the above equation. 
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( 2 ) 



Where: 


p 2 = M 2 — 1 (3) 

The terms f n (0>9n(0 in Equation 2 are singularity distributions along the body axis. The perturbation 
velocities are found by: 


8<j> 

U ~ dx 

(4) 

d<j> 

V ~ dr 

( 5 ) 

1 d4> 
W ~ r 89 

( 6 ) 

For axisymmetric flows (Reference 2) the governing P.D.E and solution are: 



(?) 

Jo - P 2 * 7 

( 8 ) 

, , r pr raw 

u(x,r) = / 

Jo v (*-0 

( 9 ) 

, , i r-* r 

v(x,r)-- 

J ° V(x- 0 -/? 2 r a 

( 10 ) 

For n = 0 the solution (Equation 2) to the non- axisymme trie problem reduces to: 


r , 

( 11 ) 

By comparison with Equation 8 it is seen that the zeroth-order solution in Equation 2 recovers the exact 
solution to the axisymmetric problem. Higher-order terms, n > 0, in Equation 2 represent non-axisymmetric 
corrections to the zeroth-order axisymmetric solution. 


2 



Determination of the Perturbation Velocities 


Returning to the non-axisymmetric problem, let: 


4>(x,r,d ) = ^2<j> n (x,r, 6) 

n=0 


Where now: 


4>n{x,r,e) = cos nd 

Jo 


x~pr | 


V 




fi r 


+ 


-^->/(x-0 2 -^r 2 

fir 




y/(x-tf-PT 


+ sin n6 


L 


x-pr i 


x-Z+J{x-tf-/Pr* 

pr 


+ 


pr 


9n (€)<*£ 




By change of variable: 


£ = x — pr cosh a 


— = — /?rsinh<7 
da 


di = 


Note that: 


r 

'x-t + y/(x- 0 2 -/? 2r2 ‘ 

n 

1 

x-Z - - p 2 r 2 

n " 


/?r 

i 

pr 

j 


The limits of integration in Equation 13 become: 


{ = 0 - ^ = cosh 1 (j^j 
£ — x — pr —+ cr — cosh” 1 (1) = 0 


Yields a new form of the perturbation potential: 



f n (x — fir cosh cr) cosh nada 


( 19 ) 


/•cosh l {xj0r) 

<Pn(x,r,e) - cos n6 / 

Jo 

y,cosh _1 (a:/0r) 

-fsin nO / g n (x — /?r cosh cr) cosh nada 

Jo 

The axial perturbation velocity may now be found by: 

oo 

u(r,r,0) = ^u n (x,r,0) 

Where: 


n = 0 


t o\ d ^ n 


( 20 ) 


( 21 ) 


cosh l (x/0r) 

fn {x “ P r cosh a) cosh nada (22) 

Q ycosh“ 1 (r/ j 0r) 

-fsi nnO— I g n (x — fir cosh cr) cosh rurda 

dx J 0 

Since the upper limit in the integrals in Equation 22 is a function of the variable of differentiation, Liebnitz 
rule* must be used. (The following applies to both of the integrals above. Only the first term will be treated 
explicitly.) 


u n (x, r, 0) = cos nO 


dx 


l 


cosh l (xfpr ) 

f n (x — fir cosh cr) cosh nada = (23) 

/•cosh" l (x/p r ) Q 

= / — [f n (x — fir cosh a) cosh na] da 

Jo ux 

g 

+/n ( x — fi r cosh cosh“ 1 (x//?r)) cosh n cosh” l (x//?r) — cosh~ l (x/ fir) 

For closed nosed bodies, / n (0) — 0, and so: 



rCOsh" l (r//3r) g 

u n (x, r, 0) = cos nO I — /„ (x — fir cosh a) cosh nada 

Jo 


+ sin n6 


l 


dx 

cosh" l (xJ0r) g 


— — g n (x — fir cosh a) cosh nada 
dx 


Recall the change of variable in Equation 14: 


j^f n (X - pr cosh <T) = (0 ^ = fn (0 


( 24 ) 


( 25 ) 


8 f iW t . . . [ k df , ... .db .. .da 
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f' n ( x — /3r cosh a) cosh ncrda 


(26) 


fcosh l (xf@r) 
D 

i>cosh _ l (xf 0r) 


u n (x , r, 0) = cos n6 

Jo 

+ sinn0 / 9n( x ~ P r cosh a ) cosh na da 

Jo 

Reversing the change of variable in Equations 14 through 18 yields. 


rx — fir 1 

u n (x t r t 6) = cos n9 - 


x-t + y/{*-tf-P 2 r 2 

pr 


P r 


n \ 


MW 




+ sin nO 


fT-Pr i 

Jo 2 


\ \x -( + \J(x-0 7 -P 2 r 2 
pr 


+ 


x-Z- \J{*-0 2 -P 2r2 

"1 


pr 

J 

J(x-i) 2 -0 7 r* 


(27) 


Similarly, the radial perturbation velocity can be found: 

oo 

v(x,r,6) = ^2v n (x,r,6) 

n= 0 


Where: 


/ n\ 

v n (x } r,0) = 


(28) 


(29) 


Q y»cosh l (x//?r) 

v n {x,r,6) = cosn0 — / /„ (* - pr cosh <r) cosh nada 

&r Jo 


g yccwh l (x/0r) 

d~rJ Q 

And by similar reasoning used in Equation 23: 


d f c 

+ sin riO-r- / g n (x - Pr cosh a) cosh nada 

dr J Q 


(30) 


^cosh _l (x//3r) g 


v n (x,r,6) = cos n6 J 5r 

fcash _l (r//?r) ^ 


f n (x - Pr cosh a) cosh nada 


+ sin nO 

Recalling the change of variable in Equation 14: 


J ( x ~~ cos ^ cos ^ ncr< * a 


fn (x - /3r cosh a) = J^/ n (0 ^ = /' (0 (~P cosh a) 


(31) 


(32) 
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/•cosh l (xfpr) 

v n (x,r, 9) = -cosn0 / ^cosher/' (x -/?r cosh <r) cosh nerefer (33) 

Jo 

/•cosh -1 (x/pr) 

— sin nO (3 cosh <t< 7 ^ (x — f3r cosh c r) cosh nada 

Jo 

Reversing the change of variable in Equations 14 through 18 yields: 



Finally, the aximuthal perturbation velocity can be found: 

oo 

w(x, r, 0) — w n (x, r, 9) 

n =0 

Where: 


And by direct differentiation: 


w n (x,r,9) 


r de 


(35) 


(36) 
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Examining the perturbation velocities for n = 0, the lowest order term in the non-axisymmetric formulation, 
it is found that: 


«o(z, 

rX-Pr 

r,6)= - 

/om 

( n 

(38) 

Jo j 

(*-{) 2 -/? 2 r 2 


t>o(x,r, 

II 

1 

1 — 

1 

-i 

(*-0/o(Otf 

(39) 

r Jo 




w 0 (x,r,6) 

= 0 

(40) 


Comparing the above with Equations 9 and 10 shows that the zeroth-order perturbation velocities recover 
the exact solution to the axisymmetric problem. Higher-order terms, n > 0, in Equations 27, 34 and 37 
represent non-axisymmetric corrections to the zeroth-order axisymmetric solution. 

Body Surface Tangent Flow Boundary Condition 

The flow along the surface of a body cannot penetrate the body, and must therefore follow the contour of 
the body. The requirement that the flow be tangent to the body, on the body surface, can be used as a 
boundary condition from which the singularity distributions may be determined. The tangent flow boundary 
condition may be written as: 


K-V5 = 0 


(41) 

Where S(x , r, 6) = 0 defines the body surface. For cylindrical coordinates in 
velocities: 

terms of the perturbation 

- - ( - \ dS _ dS xb dS 

y.V5 = ((/ + .) 57 + ^ + 7 ^ = o 


(42) 

And since u << U the boundary condition becomes: 



dS dS w dS 

dr r d9 


(43) 

Where v = v/U, etc. Defining the body surface as: 



S(x,r 9 $) = r-R{x t 0) = 0 


(44) 

Where R(x t 6) describes the body radius at a given x,6 location. It then follows that: 


OS _ dR 

dx dx 


(45) 

f-H 

II 

col^ 


(46) 

8S dR 

de “ de 


(47) 

Substituting the above into Equation 43 yields: 



dR w dR 

dx V r d6 


(48) 


Where v,w are the known velocity perturbation functions previously defined in Equations 28 and 35. The 
boundary condition is applied on the body surface, r = R(x t 0) y and so: 
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(49) 


= v(x, r, 9) - w ( x ’ r ’ 0 )R^ 0 j§Q R ( x ’ e ) 

It is convienent to introduce the following change of notation: 

u(x,r, 6 ) = 53 cos n 6 X“[f^] + 53 sin n0I" [</(,] (50) 

n=0 n= 0 

v(x,r, 6 ) = 53 cos ndX„ [/„] + 53 sin n#J” [</(,] (51) 

n=0 n=0 

oo oo 

w(x,r,e) = 53 cos ndlM - 53 sin n0X"[/„] (52) 

n=0 n=0 

Where the integral operators for each perturbation velocity are defined as: 





It is also appropriate to write R(x } 8) as a Fourier series: 


oo oo 

R(x,6) = a„(i) cos n0+ y^6 n (x)sin nO 

n=0 n=0 

(59) 


(60) 


(61) 

l^»)=|C“)-.»(l"„) 

(62) 


Substituting Equations 57, 58, and 60 through 62 into Equation 49 results in: 





(63) 


The above equation is the exact boundary condition for tangential flow on the surface of a body. The last 
term consists of products and quotients of infinite series, and presents a problem in attempting to determine 
the singularity distributions f' n ,g' n . If the body surface is very nearly axisymmetric, then dR/dQ ~ 0 and 
w « v. This limitation will be called the quasi-axisymmetric boundary condition. In this case the last 
term in Equation 63 can be dropped and the following boundary condition results: 


f cos\ / da n /dx — X“[/ n ]\ 
~i \sin/ n \db n /dx -X^\g' n ]) 


= 0 


(64) 


Note that the above boundary condition is exact for axisymmetric bodies, and is an approximation for non- 
axisymmetric bodies. In Equation 64, since cosine and sine are orthogonal functions, and since the infinite 
series are equal to zero, it necessarily follows that each term in the series must be equal to zero also. Thus 
Equation 64 reduces to the following pair of independent sets of linear equations: 


** = 0,1,2, ...oo 


TST-SW- 


n = 1,2,3, ...oo 


(65) 


(66) 


The above equations provide unique relations between the singularity distributions and the derivatives of 
the body Fourier coefficients. For known singularity distributions, a body surface is defined. For a defined 
body surface, a unique set of singularity distributions must exist. To determine the singularity distribution 
from a known body, start by expanding Equation 65. A similar analysis exists for Equation 66. 


danjx) 

dx 



_1_ 

2r 

J_ 

2r 


x -£ + y/jx -Q 2 - /? 2 r^ 

fir 

x-Z- ■\/(z-0 2 

/Jr 


(»-o/Am 

^(x-O 2 -Z ? 2 '- 2 


(67) 
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Discretization of the Singularity Distributions 

There is a relatively simple technique which can be used to determine /*(f) if da n (x)/dx is known. If /' (£) 
is assumed piecewise constant over a small interval (£*-i,£i), then /£(£;) can be pulled out of the integrals 
as follows. 


da n {x) 

dx 


£ '"«->/ .5 

1=1 

JL 1 

1=1 


i -£ + \/(*-0 2 ~^ 2 r 2 

g - £ - \A* - 0 2 - 

0r 


(x -£)(% 

\fi x - 0 2 - £ 2r2 
(»-prfe 

\ A * - o 2 - / ?2r2 


Where: 


(68) 


6 = 0 


(69) 


6 = x - /?r (70) 

Note that x,r are fixed locations on the body surface. It is convient to introduce another integral notation. 
Let: 


= /“ ± 
” It- 4., 2r 


-K* r* i 

L-1.S 

Then Equation 68 becomes: 


X -£+ ^/(g -Q 2 

0r 

' x-Z-y/jz-tf-fPr* 

0r 


(x-t) df 


J ^/( a: ~^ 2- ^ 2r2 

_ n 

(z - Q 

. \/(*-0 2 -/? 2 r 2 


da n (x) 


dx “ £/n(6) + 7 n )| fi i 

By breaking up the above summation: 


*=i 


= /;<« (;r + «•) t, +£/;<« ('»■* + 


Or finally: 


(71) 


(72) 


(73) 


(74) 


£ &&) 
ft (e,\ — * =1 


if. da n (x) 

lti-i dx 

I 

[ln + +In~) 

tk 

£fc-l 


(75) 


Thus /£(£*) can be found in terms of known integrals and /„(£fc-i), /'(fjfc_ 2 )> •••• Therefore, the entire 
singularity distribution can be determined by marching down along the body axis. The corresponding 
relation for the g f n terms is: 
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( 76 ) 


2> B «o 

_/ /<* \ _ »— i 

(/n + + 7 n") 

f. , <?6 n (x) 

ti-i dx 

3 

1! 

1 

fa + + £") 

Zk 

Zk-l 


With the singularity distributions known it is possible to calculate the pressure coefficient field on and about 
the body. A first order pressure coefficient is defined as: 

c p (x, r, 9) = -2 u(x, r, 9) ( 77 ) 

And from Equation 50: 

c p (x,r,9) = -2'^cosn91^l/^]-2j2 sinn0I nl9u] ( 78 ) 

n=0 "=0 

Since f' n ,g' n are discrete: 


c p (x,r, 9) = 


oo k y£, 1 

- 2^cosn<?^ /'(£;) / - 

n = 0 »‘= 1 

oo k y{; i 

- 2 ^2 cos n9 ^2 / 2 

n = 0 i=l 

oo k y£i 1 

- 2 sin <;(,(&) / x 

n = 0 i = l 

oo * y£i i 

- 2 sin n9 22 9n(^i) / , 

n=0 t=l J ii-\ 


'x-S+y/ix-tf-Pr*' 

dP 

Pr 

>/0 t-Pf-P 2 r 2 


n 

dP 

pr 

y/(x-P) 2 -p*r* 

'x~Z+J(x-tf-P 2 r 2 ' 

n 

dp 

pr 

xjtx-if-pr 2 

'x-S-yf(x-0 2 -P 2 r 2 ' 

n 

dp 

pr 

\J{ X — P) 2 ~ P 2 * 2 


Simplifying the above expression by the following integral definitions: 

r (> i 


jr 


* = / - 
k-x 2 


z-£ + \/( a: -0 2 -/ ?2r2 

fir 


dP 


7- = r I 

n ■«!-. i,2 

Yields the following expression for the pressure coefficient: 


*-t-^(x-Q 2 -/? 2 

pr 


J ^(x-0 2 -^ 

_ n 

y/ix-tf-Pr* 


(79) 


(80) 


(81) 


oo k 

c P (x,r,9 ) = -2^2 cos n9 ^2 f' n {Pi) 

n=0 »=1 


| f * - 2^sinn0^yJ,(£,) 

^‘ _1 n=0 i= 1 

^ + + C') 

£i-i 


Evaluation of the Integrals 

All that remains is to find closed form solutions for the integrals 7* 


£ ;-i 


r 

i ■‘n 




Ui-l 


fi 


7 U 7 U 


ii- 


ii 


(82) 


given 


in Equations 71, 72, 80, and 81. Starting with the evaluation of the integral in Equation 71, where now 
c = Pr: 
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r-f = J_/ f ’ 

n l fi _, 2rc n Je 


^ ^ \/(x — £) — c 2 


Introducing the first change of variable: 


<7 = x — £ 


Yields: 


,: 1L = '2^/,J' T + 


n <7d<7 


Introducing the second change of variable: 


p — vV 2 - c 2 


j “ c » 

a <7 — dp 

a 


Yields: 




Introducing the third change of variable, and utilizing one of the indefinite integrals found in the CRC 
Standard Math Tables* : 


A = tan 1 - 


Yields: 


7„ + — — C / [ctan A -f csec A] n sec 2 AdA 

f<-i 2rc n y A ._ x 

v + f 1 -h sin A 1 n dA 

n *.-i 2r/ x ._ 1 [ cos A J cos 2 A 


Introducing the fourth, and final, change of variable: 


1 + sin A 


d\ = cos A — 


«+ ^ r ^ A 

n *-i 2rJ rii ^ cos A ’ 


J f {z> ^ z 2 + fl2 ) dx = a J f (a tan u, a sec u) sec 2 udu, u = tan 1 - , a > 0 
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From Equation 92, and utilizing simple trigonometric identities, it can be shown that: 


cos A = 


T] 2 + 1 


Substituting Equation 95 into Equation 94, and recalling that c = pr, yields a simple expression which can 
be directly integrated in closed form: 


For n = 0: 


For n = 1: 


For n > 2: 


nt =-§ f ■r’tf+O 

^ Jqi-i 


el" T 


" If.-. 4 l »J„_, 


,: X.r +x '"' 


r + | e ’ = -(L\J—T 1 n + 1 + -!—T 1 n - 1 Y' (99) 

n \u-x 4[n+l' n-1 J„._, 

Reversing the change of variables in Equations 92, 89, 86, and 84 yields the limits of integration in Equations 
97, 98, and 99: 


■-ft-i + xAz-fc- 1) 2 -/? 2 ” 2 


tt-i = 




Continuing on with the evaluation of the integral in Equation 72: 


” «<-. 2"" 4., I V ' J 


l£»-i ^ rC L J V( x “0 *“ c2 

Using the same change of variable in Equations 84, 86, and 89 yields the following relations. 


n ada 


J® ^ = f [c sec A — c tan A] n sec 2 Ad A 

C«-i 2rc n Jxi^x 

c f Ai ri—sinAl n dA 
7 " Ci-1 ~ 2r A.., [ cos A J cos 2 A 


f Ai 1— sinAl n dA 
> cos A cos 2 A 


Introducing a new fourth change of variable: 


1 — sin A 
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Yields: 


d A = — cos A — 
P 


(108) 


II 


From Equation 107, it can be shown that: 




.n — 1 


COS A 


dp 


(109) 


cos A = 


2p 

M 2 + 1 


( 110 ) 


Substituting Equation 110 into Equation 109 yields a simple expression which can be directly integrated in 
closed form: 


For n — 0: 


For n = 1: 



§ r /.-mp’+ihk 

q J m-i 


(in) 


r' 

1 n 


Zi 

Zi-l 



(112) 


For n > 2: 



£ 

4 


lp 2 + ln/i 


(113) 




* 

* „n+i . A_ /i”- 1 

1^ 

1 

7 

. ! P T 7 P 

n -f 1 n — 1 


A*«-l 


(114) 


Reversing the change of variables in Equations 107, 89, 86, and 84 yields the limits of integration in Equations 
112, 113, and 114: 


Pi-i 


_ X-Si-l - ^/( z - 6 - i) 2 -/? 2 t 


Pr 


x-£i- J{x- iif - P 2 

W " Tr 

Proceding with the evaluation of the integral in Equation 80: 


if 


L. - =? jC. [— «+Vc— o’- 




\/(*-0 2 -c 2 

Using the same change of variable in Equations 84, 86, 89, and 92 yields the following: 

da 


rr 


In 


L. = _ S 


V<r 2 - c 2 

_fp == 

s/p 1 + c 1 


(115) 

(116) 

(117) 

(118) 
(119) 
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(120) 


= f [c tan A + csec A] n sec XdX 

U«-i 2c n Jxi-i 

U 1 f Xi I* 1 + sin A] 

fi-1 2 1 L 


cos A 


n dX 
cos A 


/ u 

1 n 


• _ i r 

2 J ^ 

‘-i * Jfji-i 




For n = 0: 


For n = 1: 


For n > 2: 


7 U 

■* n 


/r 


e. 

Ci-1 

\Zi 

le,-i 


\\nrj 


J t 


= - o'? 


1«J. 




if 



' 1 

1 

II 

1 

2n^ . 


The limits of integration for Equations 123, 124, and 125 are the same as those given in Equations 
101. Finally, the evaluation of the last integral in Equation 81: 








\/(x-£f-c 2 

Using the same change of variable in Equations 84, 86, 89, and 107 yields the following. 


I u ~ 

r~ 

* n 


dp 

\Jp 2 + 1 


if 


= — f [csec A — ctan Al" 

<-> 2c" A.., 

< _ _i ( Xi f l- si 

c-i : =\r»'-' 

k.-i z Jm-\ 


sec XdX 
in dX 


dfi 


For n = 0: 


For n = 1: 


/ u " 




l£»-i 




if 


£•-1 


= 






^i-i 


(121) 

( 122 ) 

(123) 

(124) 

(125) 
100 and 

(126) 

(127) 

(128) 

(129) 

(130) 

(131) 

(132) 

(133) 
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For n > 2: 


/ u 
•* n 


1 

1 "1 

11 

7 

[to*] 






(134) 


The limits of integration for Equations 132, 133, and 134 are the same as those given in Equations 115 
and 116. This completes the evaluation of all the integrals necessary to formulate the problem solution 
numerically. 


Summary 

An analysis technique has been developed employing supersonic linearized small perturbation theory which 
is not restricted to slender bodies of circular cross section. The method developed here can be employed in 
the analysis of wave drag for slender, arbitrary cross section bodies. The method permits the calculation of 
the first-order pressure coefficient on or about the body. The pressure coefficient can be integrated over the 
surface of the body to compute the wave drag resulting from the body. Also, the pressure coefficient field 
can be calculated about the body, and this pressure field can then be used in other existing analyses and 
computer programs to determine the aerodynamic interference drag for complex aircraft configurations. 
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Figure 1; Cylindrical Coordinate System 
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